Relatedness and genomic inbreeding in a sample of Timor ponies.

Timor ponies (TP) were first shipped to Australia in the early 1800s and were highly valued as transport and pack animals, which resulted in TPs contributing to the development of Australian horse breeds. Today, while the exact number of TPs in Australia is currently unknown, there has been recent interest in establishing a domestic breeding program for Australian TPs. The aim of this study was to evaluate the relatedness of a sample of TPs, as well as provide estimates of genomic inbreeding levels to better inform the feasibility of using these animals as founders for a domestic breeding program. Hair samples from each horse were genotyped using the Illumina 80K Infinium Equine genotyping array and data were analysed using PLINK v1.90b7, KING 2.3.2 and R v4.3.1. The results illustrate that there are distantly related and minimally inbred horses within the sampled TPs. Lengths of the ROH segments also indicated that recent inbreeding events are likely to only have occurred in a third of the horses. Overall, these results are promising for the success of a domestic TP breeding program; however, considering the low number of domestic TPs known to reside in Australia, there would certainly still be substantial benefits to incorporating additional TPs either directly from Timor or from areas in Australia that are believed to contain wild descendants of TPs.

Epichaperome inhibition targets TP53-mutant AML and AML stem/progenitor cells.

TP 53-mutant acute myeloid leukemia (AML) remains the ultimate therapeutic challenge. Epichaperomes, formed in malignant cells, consist of heat shock protein 90 (HSP90) and associated proteins that support the maturation, activity, and stability of oncogenic kinases and transcription factors including mutant p53. High-throughput drug screening identified HSP90 inhibitors as top hits in isogenic TP53-wild-type (WT) and -mutant AML cells. We detected epichaperomes in AML cells and stem/progenitor cells with TP53 mutations but not in healthy bone marrow (BM) cells. Hence, we investigated the therapeutic potential of specifically targeting epichaperomes with PU-H71 in TP53-mutant AML based on its preferred binding to HSP90 within epichaperomes. PU-H71 effectively suppressed cell intrinsic stress responses and killed AML cells, primarily by inducing apoptosis; targeted TP53-mutant stem/progenitor cells; and prolonged survival of TP53-mutant AML xenograft and patient-derived xenograft models, but it had minimal effects on healthy human BM CD34+ cells or on murine hematopoiesis. PU-H71 decreased MCL-1 and multiple signal proteins, increased proapoptotic Bcl-2-like protein 11 levels, and synergized with BCL-2 inhibitor venetoclax in TP53-mutant AML. Notably, PU-H71 effectively killed TP53-WT and -mutant cells in isogenic TP53-WT/TP53-R248W Molm13 cell mixtures, whereas MDM2 or BCL-2 inhibition only reduced TP53-WT but favored the outgrowth of TP53-mutant cells. Venetoclax enhanced the killing of both TP53-WT and -mutant cells by PU-H71 in a xenograft model. Our data suggest that epichaperome function is essential for TP53-mutant AML growth and survival and that its inhibition targets mutant AML and stem/progenitor cells, enhances venetoclax activity, and prevents the outgrowth of venetoclax-resistant TP53-mutant AML clones. These concepts warrant clinical evaluation.

Reconstruction of the Major Maternal and Paternal Lineages in the Feral New Zealand Kaimanawa Horses.

New Zealand has the fourth largest feral horse population in the world. The Kaimanawas (KHs) are feral horses descended from various domestic horse breeds released into the Kaimanawa ranges in the 19th and 20th centuries. Over time, the population size has fluctuated dramatically due to hunting, large-scale farming and forestry. Currently, the herd is managed by an annual round-up, limiting the number to 300 individuals to protect the native ecosystem. Here, we genotyped 96 KHs for uniparental markers (mitochondrial DNA, Y-chromosome) and assessed their genetic similarity with respect to other domestic horses. We show that at least six maternal and six paternal lineages contributed unequally to the KH gene pool, and today's KH population possibly represents two sub-populations. Our results indicate that three horse breeds, namely Welsh ponies, Thoroughbreds and Arabian horses had a major influence in the genetic-makeup of the extant KH population. We show that mitochondrial genetic diversity in KHs (π = 0.00687 ± 0.00355) is closer to that of the Sable Island horses (π = 0.0034 ± 0.00301), and less than other feral horse populations around the world. Our current findings, combined with ongoing genomic research, will provide insight into the population-specific genetic variation and inbreeding among KHs. This will largely advance equine research and improve the management of future breeding programs of these treasured New Zealand horse.

Y Chromosome Haplotypes Enlighten Origin, Influence, and Breeding History of North African Barb Horses.

In horses, demographic patterns are complex due to historical migrations and eventful breeding histories. Particularly puzzling is the ancestry of the North African horse, a founding horse breed, shaped by numerous influences throughout history. A genetic marker particularly suitable to investigate the paternal demographic history of populations is the non-recombining male-specific region of the Y chromosome (MSY). Using a recently established horse MSY haplotype (HT) topology and KASP™ genotyping, we illustrate MSY HT spectra of 119 Barb and Arab-Barb males, collected from the Maghreb region and European subpopulations. All detected HTs belonged to the Crown haplogroup, and the broad MSY spectrum reflects the wide variety of influential stallions throughout the breed's history. Distinct HTs and regional disparities were characterized and a remarkable number of early introduced lineages were observed. The data indicate recent refinement with Thoroughbred and Arabian patrilines, while 57% of the dataset supports historical migrations between North Africa and the Iberian Peninsula. In the Barb horse, we detected the HT linked to Godolphin Arabian, one of the Thoroughbred founders. Hence, we shed new light on the question of the ancestry of one Thoroughbred patriline. We show the strength of the horse Y chromosome as a genealogical tool, enlighten recent paternal history of North African horses, and set the foundation for future studies on the breed and the formation of conservation breeding programs.

Copy number variation of horse Y chromosome genes in normal equine populations and in horses with abnormal sex development and subfertility: relationship of copy number variations with Y haplogroups.

Structural rearrangements like copy number variations in the male-specific Y chromosome have been associated with male fertility phenotypes in human and mouse but have been sparsely studied in other mammalian species. Here, we designed digital droplet PCR assays for 7 horse male-specific Y chromosome multicopy genes and SRY and evaluated their absolute copy numbers in 209 normal male horses of 22 breeds, 73 XY horses with disorders of sex development and/or infertility, 5 Przewalski's horses and 2 kulans. This established baseline copy number for these genes in horses. The TSPY gene showed the highest copy number and was the most copy number variable between individuals and breeds. SRY was a single-copy gene in most horses but had 2-3 copies in some indigenous breeds. Since SRY is flanked by 2 copies of RBMY, their copy number variations were interrelated and may lead to SRY-negative XY disorders of sex development. The Przewalski's horse and kulan had 1 copy of SRY and RBMY. TSPY and ETSTY2 showed significant copy number variations between cryptorchid and normal males (P < 0.05). No significant copy number variations were observed in subfertile/infertile males. Notably, copy number of TSPY and ETSTY5 differed between successive male generations and between cloned horses, indicating germline and somatic mechanisms for copy number variations. We observed no correlation between male-specific Y chromosome gene copy number variations and male-specific Y chromosome haplotypes. We conclude that the ampliconic male-specific Y chromosome reference assembly has deficiencies and further studies with an improved male-specific Y chromosome assembly are needed to determine selective constraints over horse male-specific Y chromosome gene copy number and their relation to stallion reproduction and male biology.

Y-Chromosomal Insights into Breeding History and Sire Line Genealogies of Arabian Horses.

The Y chromosome is a valuable genetic marker for studying the origin and influence of paternal lineages in populations. In this study, we conducted Y-chromosomal lineage-tracing in Arabian horses. First, we resolved a Y haplotype phylogeny based on the next generation sequencing data of 157 males from several breeds. Y-chromosomal haplotypes specific for Arabian horses were inferred by genotyping a collection of 145 males representing most Arabian sire lines that are active around the globe. These lines formed three discrete haplogroups, and the same haplogroups were detected in Arabian populations native to the Middle East. The Arabian haplotypes were clearly distinct from the ones detected in Akhal Tekes, Turkoman horses, and the progeny of two Thoroughbred foundation sires. However, a haplotype introduced into the English Thoroughbred by the stallion Byerley Turk (1680), was shared among Arabians, Turkomans, and Akhal Tekes, which opens a discussion about the historic connections between Oriental horse types. Furthermore, we genetically traced Arabian sire line breeding in the Western World over the past 200 years. This confirmed a strong selection for relatively few male lineages and uncovered incongruences to written pedigree records. Overall, we demonstrate how fine-scaled Y-analysis contributes to a better understanding of the historical development of horse breeds.

Investigations on the occurrence of a muscular disorder in Austrian slaughter pigs.

BACKGROUND: In slaughterhouse, veterinarians responsible for meat inspection are often confronted with abnormalities in carcasses, not only in pigs but in all livestock species. In 2017, a veterinarian responsible for meat inspection in a slaughterhouse in Styria, Austria, observed gluteal muscles infiltrated by fat and muscle tissue obviously being replaced by fat in two different slaughter pigs. These muscles were sent for further diagnostics to the University Clinic for Swine. RESULTS: The two muscle samples were investigated histopathologically and diagnosed with fatty muscular dystrophy. The results of routine histopathology were confirmed by dystrophin-specific immunohistochemistry. Sex of the two affected animals was determined retrospectively using a PCR-based protocol and resulted in one male and one female pig. A survey to determine the prevalence of fatty muscular disorders of pork revealed that this phenomenon gets frequently observed in Styria, but also occurs in Upper Austria and Lower Austria. Mostly gluteal and lumbal muscles were affected and approximately 20-40% of the affected muscles were replaced by fat. CONCLUSIONS: Fatty muscular dystrophy or muscular steatosis, as it was sometimes called in early literature, seems not to be an uncommon and rare event and is known to have several different causes. As it was detected in both sexes, our observations are different to the described case in Japan, where only one male individual was affected. To avoid further increase of such cases (fatty muscular dystrophy), it would be useful to clarify the cause. First, whether the cause is environmental or genetic, and in case it is genetic it would be key to disentangle the underlying genomic architecture. Having causal variants described-one could think about integrating this information (depending on the mode of inheritance and the number of loci involved) in the breeding program of pigs. Furthermore, the proportion of non-Austrian pig genetics used for commercial pig production in Austria should be reviewed in order to be able to make reliable statements about the spread of the disease not only in Austrian pig breeds, but also in pig breeds worldwide.

Distribution of the Warmblood Fragile Foal Syndrome Type 1 Mutation (PLOD1 c.2032G>A) in Different Horse Breeds from Europe and the United States.

Warmblood fragile foal syndrome (WFFS) is an autosomal recessive disorder caused by a single nucleotide variant in the procollagen-lysine-2-oxoglutarate-5-dioxygenase 1 gene (PLOD1:c.2032G>A, p.Gly678Arg). Homozygosity for the PLOD1 variant causes an Ehler-Danlos-like syndrome, which has to date only been reported in warmblood breeds but the WFFS allele has been also detected in the Thoroughbred. To investigate the breed distribution of the WFFS allele, 4081 horses belonging to 38 different breeds were screened. In total, 4.9% of the horses representing 21 breeds carried the WFFS allele. The affected breeds were mainly warmbloods, with carrier frequency as high as 17% in the Hanoverian and Danish Warmblood. The WFFS allele was not detected in most non-warmblood breeds. Exceptions include WFFS carriers in the Thoroughbred (17/716), Haflinger (2/48), American Sport Pony (1/12), and Knabstrupper (3/46). The origin of the WFFS allele remains unknown. The Arabian breed and specifically the stallion Bairactar Or. Ar. (1813), whose offspring were reported to have a similar phenotype in the 19th century, were hypothesized as the origin. DNA from a museum sample of Bairactar Or. Ar. showed that he did not carry the mutated allele. This result, together with the genotypes of 302 Arabians, all homozygous for the reference allele, does not support an Arabian origin of the WFFS allele. Our extensive survey shows the WFFS allele to be of moderate frequency and concern in warmbloods and also in breeds where it may not be expected.

Genome Diversity and the Origin of the Arabian Horse.

The Arabian horse, one of the world's oldest breeds of any domesticated animal, is characterized by natural beauty, graceful movement, athletic endurance, and, as a result of its development in the arid Middle East, the ability to thrive in a hot, dry environment. Here we studied 378 Arabian horses from 12 countries using equine single nucleotide polymorphism (SNP) arrays and whole-genome re-sequencing to examine hypotheses about genomic diversity, population structure, and the relationship of the Arabian to other horse breeds. We identified a high degree of genetic variation and complex ancestry in Arabian horses from the Middle East region. Also, contrary to popular belief, we could detect no significant genomic contribution of the Arabian breed to the Thoroughbred racehorse, including Y chromosome ancestry. However, we found strong evidence for recent interbreeding of Thoroughbreds with Arabians used for flat-racing competitions. Genetic signatures suggestive of selective sweeps across the Arabian breed contain candidate genes for combating oxidative damage during exercise, and within the "Straight Egyptian" subgroup, for facial morphology. Overall, our data support an origin of the Arabian horse in the Middle East, no evidence for reduced global genetic diversity across the breed, and unique genetic adaptations for both physiology and conformation.

A Genome-Wide Association Analysis in Noriker Horses Identifies a SNP Associated With Roan Coat Color.

The roan coat color in horses is characterized by dispersed white hair and dark points. This phenotype segregates in a broad range of horse breeds, while the underlying genetic background is still unknown. Previous studies mapped the roan locus to the KIT gene on equine chromosome 3 (ECA3). However, this association could not be validated across different horse breeds. Performing a genome-wide association analysis (GWAS) in Noriker horses, we identified a single nucleotide polymorphism (SNP) (ECA3:g.79,543.439 A > G) in the intron 17 of the KIT gene. The G -allele of the top associated SNP was present in other roan horses, namely Quarter Horse, Murgese, Slovenian, and Belgian draught horse, while it was absent in a panel of 15 breeds, including 657 non-roan horses. In further 379 gray Lipizzan horses, eight animals exhibited a heterozygous genotype (A/G). Comparative whole-genome sequence analysis of the KIT region revealed two deletions in the downstream region (ECA3:79,533,217_79,533,224delTCGTCTTC; ECA3:79,533,282_79,533,285delTTCT) and a 3 bp deletion combined with 17 bp insertion in intron 20 of KIT (ECA3:79,588,128_79,588,130delinsTTATCTCTATAGTAGTT). Within the Noriker sample, these loci were in complete linkage disequilibrium (LD) with the identified top SNP. Based upon pedigree information and historical records, we were able to trace back the genetic origin of roan coat color to a baroque gene pool. Furthermore, our data suggest allelic heterogeneity and the existence of additional roan alleles in ponies and breeds related to the English Thoroughbred. In order to study the roan phenotype segregating in those breeds, further association and verification studies are required.

Population Genetic Analysis of the Estonian Native Horse Suggests Diverse and Distinct Genetics, Ancient Origin and Contribution from Unique Patrilines.

The Estonian Native Horse (ENH) is a medium-size pony found mainly in the western islands of Estonia and is well-adapted to the harsh northern climate and poor pastures. The ancestry of the ENH is debated, including alleged claims about direct descendance from the extinct Tarpan. Here we conducted a detailed analysis of the genetic makeup and relationships of the ENH based on the genotypes of 15 autosomal short tandem repeats (STRs), 18 Y chromosomal single nucleotide polymorphisms (SNPs), mitochondrial D-loop sequence and lateral gait allele in DMRT3. The study encompassed 2890 horses of 61 breeds, including 33 ENHs. We show that the expected and observed genetic diversities of the ENH are among the highest within 52 global breeds, and the highest among 8 related Northern European ponies. The genetically closest breeds to the ENH are the Finn Horse, and the geographically more distant primitive Hucul and Konik. ENH matrilines are diverse and relate to draught and Pontic-Caspian breeds. ENH patrilines relate to draught breeds, and to a unique haplogroup not described before. None of the 33 ENHs carried the "gait" mutation, but the mutation was found in 2 Huculs. The study demonstrates that the ENH is a genetically distinct and diverse breed of ancient origin with no notable pressure of selective breeding.

A First Y-Chromosomal Haplotype Network to Investigate Male-Driven Population Dynamics in Domestic and Wild Bactrian Camels.

Polymorphic markers on the male-specific part of the Y chromosome (MSY) provide useful information for tracking male genealogies. While maternal lineages are well studied in Old World camelids using mitochondrial DNA, the lack of a Y-chromosomal reference sequence hampers the analysis of male-driven demographics. Recently, a shotgun assembly of the horse MSY was generated based on short read next generation sequencing data. The haplotype network resulting from single copy MSY variants using the assembly as a reference revealed sufficient resolution to trace individual male lines in this species. In a similar approach we generated a 3.8 Mbp sized assembly of the MSY of Camelus bactrianus. The camel MSY assembly was used as a reference for variant calling using short read data from eight Old World camelid individuals. Based on 596 single nucleotide variants we revealed a Y-phylogenetic network with seven haplotypes. Wild and domestic Bactrian camels were clearly separated into two different haplogroups with an estimated divergence time of 26,999 ± 2,268 years. Unexpectedly, one wild camel clustered into the domestic Bactrian camels' haplogroup. The observation of a domestic paternal lineage within the wild camel population is concerning in view of the importance to conserve the genetic integrity of these highly endangered species in their natural habitat.

The horse Y chromosome as an informative marker for tracing sire lines.

Analysis of the Y chromosome is the best-established way to reconstruct paternal family history in humans. Here, we applied fine-scaled Y-chromosomal haplotyping in horses with biallelic markers and demonstrate the potential of our approach to address the ancestry of sire lines. We de novo assembled a draft reference of the male-specific region of the Y chromosome from Illumina short reads and then screened 5.8 million basepairs for variants in 130 specimens from intensively selected and rural breeds and nine Przewalski's horses. Among domestic horses we confirmed the predominance of a young'crown haplogroup' in Central European and North American breeds. Within the crown, we distinguished 58 haplotypes based on 211 variants, forming three major haplogroups. In addition to two previously characterised haplogroups, one observed in Arabian/Coldblooded and the other in Turkoman/Thoroughbred horses, we uncovered a third haplogroup containing Iberian lines and a North African Barb Horse. In a genealogical showcase, we distinguished the patrilines of the three English Thoroughbred founder stallions and resolved a historic controversy over the parentage of the horse 'Galopin', born in 1872. We observed two nearly instantaneous radiations in the history of Central and Northern European Y-chromosomal lineages that both occurred after domestication 5,500 years ago.

Ancient genomes revisit the ancestry of domestic and Przewalski's horses.

The Eneolithic Botai culture of the Central Asian steppes provides the earliest archaeological evidence for horse husbandry, ~5500 years ago, but the exact nature of early horse domestication remains controversial. We generated 42 ancient-horse genomes, including 20 from Botai. Compared to 46 published ancient- and modern-horse genomes, our data indicate that Przewalski's horses are the feral descendants of horses herded at Botai and not truly wild horses. All domestic horses dated from ~4000 years ago to present only show ~2.7% of Botai-related ancestry. This indicates that a massive genomic turnover underpins the expansion of the horse stock that gave rise to modern domesticates, which coincides with large-scale human population expansions during the Early Bronze Age.

Exchange of genetic information between therian X and Y chromosome gametologs in old evolutionary strata.

Therian X and Y sex chromosomes arose from a pair of autosomes. Y chromosomes consist of a pseudoautosomal region that crosses over with the X chromosome and a male-specific Y-chromosomal region that does not. The X chromosome can be structured into "evolutionary strata". Divergence of X-chromosomal genes from their gametologs is similar within a stratum, but differs among strata, likely caused by a different onset of suppression of crossing over between gametologs. After stratum formation, exchange of information between gametologs has long been believed absent; however, recent studies have shown limited exchange, likely through gene conversion. Herein we investigate exchange of genetic information between gametologs in old strata that formed before the split of Laurasiatheria (cattle) from Euarchontoglires (primates and rodents) with a new phylogenetic approach. A prerequisite for our test is an overall preradiative topology, that is, all X-chromosomal gametologs are more similar among themselves than to Y-chromosomal sequences. Screening multiple sequence alignments of the coding sequences of genes from cattle, mice, and humans identified four genes, DDX3X/Y,RBMX/Y,USP9X/Y, and UTX/Y, exhibiting a preradiation topology. Applying our test, we detected exchange of genetic information between all four X and Y gametologs after stratum formation.

Developing a 670k genotyping array to tag ~2M SNPs across 24 horse breeds.

BACKGROUND: To date, genome-scale analyses in the domestic horse have been limited by suboptimal single nucleotide polymorphism (SNP) density and uneven genomic coverage of the current SNP genotyping arrays. The recent availability of whole genome sequences has created the opportunity to develop a next generation, high-density equine SNP array. RESULTS: Using whole genome sequence from 153 individuals representing 24 distinct breeds collated by the equine genomics community, we cataloged over 23 million de novo discovered genetic variants. Leveraging genotype data from individuals with both whole genome sequence, and genotypes from lower-density, legacy SNP arrays, a subset of ~5 million high-quality, high-density array candidate SNPs were selected based on breed representation and uniform spacing across the genome. Considering probe design recommendations from a commercial vendor (Affymetrix, now Thermo Fisher Scientific) a set of ~2 million SNPs were selected for a next-generation high-density SNP chip (MNEc2M). Genotype data were generated using the MNEc2M array from a cohort of 332 horses from 20 breeds and a lower-density array, consisting of ~670 thousand SNPs (MNEc670k), was designed for genotype imputation. CONCLUSIONS: Here, we document the steps taken to design both the MNEc2M and MNEc670k arrays, report genomic and technical properties of these genotyping platforms, and demonstrate the imputation capabilities of these tools for the domestic horse.

Y Chromosome Uncovers the Recent Oriental Origin of Modern Stallions.

The Y chromosome directly reflects male genealogies, but the extremely low Y chromosome sequence diversity in horses has prevented the reconstruction of stallion genealogies [1, 2]. Here, we resolve the first Y chromosome genealogy of modern horses by screening 1.46 Mb of the male-specific region of the Y chromosome (MSY) in 52 horses from 21 breeds. Based on highly accurate pedigree data, we estimated the de novo mutation rate of the horse MSY and showed that various modern horse Y chromosome lineages split much later than the domestication of the species. Apart from few private northern European haplotypes, all modern horse breeds clustered together in a roughly 700-year-old haplogroup that was transmitted to Europe by the import of Oriental stallions. The Oriental horse group consisted of two major subclades: the Original Arabian lineage and the Turkoman horse lineage. We show that the English Thoroughbred MSY was derived from the Turkoman lineage and that English Thoroughbred sires are largely responsible for the predominance of this haplotype in modern horses.

A missense mutation in damage-specific DNA binding protein 2 is a genetic risk factor for limbal squamous cell carcinoma in horses.

Squamous cell carcinoma (SCC) is the most common cancer of the equine eye, frequently originating at the limbus, with the potential to invade the cornea, cause visual impairment, and result in loss of the eye. Several breeds of horses have a high occurrence of limbal SCC implicating a genetic basis for limbal SCC predisposition. Pedigree analysis in the Haflinger breed supports a simple recessive mode of inheritance and a genome-wide association study (N = 23) identified a 1.5 Mb locus on ECA12 significantly associated with limbal SCC (Pcorrected = 0.04). Sequencing the most physiologically relevant gene from this locus, damage specific DNA binding protein 2 (DDB2), identified a missense mutation (c.1013 C > T p.Thr338Met) that was strongly associated with limbal SCC (P = 3.41 × 10-10 ). Genotyping 42 polymorphisms narrowed the ECA12 candidate interval to 483 kb but did not identify another variant that was more strongly associated. DDB2 binds to ultraviolet light damaged DNA and recruits other proteins to perform global genome nucleotide excision repair. Computational modeling predicts this mutation to be deleterious by altering conformation of the ß loop involved in photolesion recognition. This DDB2 variant was also detected in two other closely related breeds with reported cases of ocular SCC, the Belgian and the Percheron, suggesting it may also be a SCC risk factor in these breeds. Furthermore, in humans xeroderma pigmentosum complementation group E, a disease characterized by sun sensitivity and increased risk of cutaneous SCC and melanomas, is explained by mutations in DDB2. Cross-species comparison remains to be further evaluated.

High-throughput mRNA and miRNA profiling of epithelial-mesenchymal transition in MDCK cells.

BACKGROUND: Epithelial-mesenchymal transition (EMT) is an important process in embryonic development, especially during gastrulation and organ formation. Furthermore EMT is widely observed in pathological conditions, e.g., fibrosis, tumor progression and metastasis. Madin-Darby Canine Kidney (MDCK) cells are widely used for studies of EMT and epithelial plasticity. MDCK cells show an epithelial phenotype, while oncogenic Ras-transformed MDCK (MDCK-Ras) cells undergo EMT and show a mesenchymal phenotype. METHODS: RNA-Seq and miRNA-Seq analyses were performed on MDCK and MDCK-Ras cells. Data were validated by qRT-PCR. Gene signature analyses were carried out to identify pathways and gene ontology terms. For selected miRNAs target prediction was performed. RESULTS: With RNA-Seq, mRNAs of approximately half of the genes known for dog were detected. These were screened for differential regulation during Ras-induced EMT. We went further and performed gene signature analyses and found Gene Ontology (GO) terms and pathways important for epithelial polarity and implicated in EMT. Among the identified pathways, TGFß1 emerged as a central signaling factor in many EMT related pathways and biological processes. With miRNA-Seq, approximately half of the known canine miRNAs were found expressed in MDCK and MDCK-Ras cells. Furthermore, among differentially expressed miRNAs, miRNAs that are known to be important regulators of EMT were detected and new candidates were predicted. New dog miRNAs were discovered after aligning our reads to that of other species in miRBase. Importantly, we could identify 25 completely novel miRNAs with a stable hairpin structure. Two of these novel miRNAs were differentially expressed. We validated the two novel miRNAs with the highest read counts by RT-qPCR. Target prediction of a particular novel miRNA highly expressed in mesenchymal MDCK-Ras cells revealed that it targets components of epithelial cell junctional complexes. Combining target prediction for the most upregulated miRNAs and validation of the targets in MDCK-Ras cells with pathway analysis allowed us to identify two novel pathways, e.g., JAK/STAT signaling and pancreatic cancer pathways. These pathways could not be detected solely by gene set enrichment analyses of RNA-Seq data. CONCLUSION: With deep sequencing data of mRNAs and miRNAs of MDCK cells and of Ras-induced EMT in MDCK cells, differentially regulated mRNAs and miRNAs are identified. Many of the identified genes are within pathways known to be involved in EMT. Novel differentially upregulated genes in MDCK cells are interferon stimulated genes and genes involved in Slit and Netrin signaling. New pathways not yet linked to these processes were identified. A central pathway in Ras induced EMT is TGFß signaling, which leads to differential regulation of many target genes, including miRNAs. With miRNA-Seq we identified miRNAs involved in either epithelial cell biology or EMT. Finally, we describe completely novel miRNAs and their target genes.